Temperature-driven coordination of circadian transcriptional regulation.

The circadian clock is an evolutionarily-conserved molecular oscillator that enables species to anticipate rhythmic changes in their environment. At a molecular level, the core clock genes induce circadian oscillations in thousands of genes in a tissue-specific manner, orchestrating myriad biological processes. While previous studies have investigated how the core clock circuit responds to environmental perturbations such as temperature, the downstream effects of such perturbations on circadian regulation remain poorly understood. By analyzing bulk-RNA sequencing of Drosophila fat bodies harvested from flies subjected to different environmental conditions, we demonstrate a highly condition-specific circadian transcriptome: genes are cycling in a temperature-specific manner, and the distributions of their phases also differ between the two conditions. Further employing a reference-based gene regulatory network (Reactome), we find evidence of increased gene-gene coordination at low temperatures and synchronization of rhythmic genes that are network neighbors. We report that the phase differences between cycling genes increase as a function of geodesic distance in the low temperature condition, suggesting increased coordination of cycling on the gene regulatory network. Our results suggest a potential mechanism whereby the circadian clock mediates the fly's response to seasonal changes in temperature.

Detecting Rhythmic Gene Expression in Single Cell Transcriptomics.

An autonomous, environmentally-synchronizable circadian rhythm is a ubiquitous feature of life on Earth. In multicellular organisms, this rhythm is generated by a transcription-translation feedback loop present in nearly every cell that drives daily expression of thousands of genes in a tissue-dependent manner. Identifying the genes that are under circadian control can elucidate the mechanisms by which physiological processes are coordinated in multicellular organisms. Today, transcriptomic profiling at the single-cell level provides an unprecedented opportunity to understand the function of cell-level clocks. However, while many cycling detection algorithms have been developed to identify genes under circadian control in bulk transcriptomic data, it is not known how best to adapt these algorithms to single-cell RNAseq data. Here, we benchmark commonly used circadian detection methods on their reliability and efficiency when applied to single cell RNAseq data. Our results provide guidance on adapting existing cycling detection methods to the single-cell domain, and elucidate opportunities for more robust and efficient rhythm detection in single-cell data. We also propose a subsampling procedure combined with harmonic regression as an efficient, reliable strategy to detect circadian genes in the single-cell setting.

Temperature-driven coordination of circadian transcriptome regulation.

The circadian rhythm is an evolutionarily-conserved molecular oscillator that enables species to anticipate rhythmic changes in their environment. At a molecular level, the core clock genes induce a circadian oscillation in thousands of genes in a tissue-specific manner, orchestrating myriad biological processes. While studies have investigated how the core clock circuit responds to environmental perturbations such as temperature, the downstream effects of such perturbations on circadian regulation remain poorly understood. By analyzing bulk-RNA sequencing of Drosophila fat bodies harvested from flies subjected to different environmental conditions, we demonstrate a highly condition-specific circadian transcriptome. Further employing a reference-based gene regulatory network (Reactome), we find evidence of increased gene-gene coordination at low temperatures and synchronization of rhythmic genes that are network neighbors. Our results point to the mechanisms by which the circadian clock mediates the fly's response to seasonal changes in temperature.

A dynamical systems treatment of transcriptomic trajectories in hematopoiesis.

Inspired by Waddington's illustration of an epigenetic landscape, cell-fate transitions have been envisioned as bifurcating dynamical systems, wherein exogenous signaling dynamics couple to the enormously complex signaling and transcriptional machinery of a cell to elicit qualitative transitions in its collective state. Single-cell RNA sequencing (scRNA-seq), which measures the distributions of possible transcriptional states in large populations of differentiating cells, provides an alternate view, in which development is marked by the variations of a myriad of genes. Here, we present a mathematical formalism for rigorously evaluating, from a dynamical systems perspective, whether scRNA-seq trajectories display statistical signatures consistent with bifurcations and, as a case study, pinpoint regions of multistability along the neutrophil branch of hematopoeitic differentiation. Additionally, we leverage the geometric features of linear instability to identify the low-dimensional phase plane in gene expression space within which the multistability unfolds, highlighting novel genetic players that are crucial for neutrophil differentiation. Broadly, we show that a dynamical systems treatment of scRNA-seq data provides mechanistic insights into the high-dimensional processes of cellular differentiation, taking a step toward systematic construction of mathematical models for transcriptomic dynamics.

LPS-Induced Systemic Inflammation Caused mPOA-FSH/LH Disturbance and Impaired Testicular Function.

Male reproductive function is key to the continuation of species and is under sophisticated regulation, challenged by various stressors including inflammation. In the lipopolysaccharide (LPS) intraperitoneal injection-induced acute systemic inflammation, male fecundity was compromised with decreased testosterone level, damaged spermatogenesis, and downregulations of testicular gene expression levels involved in steroidogenesis regulation and blood-testis barrier. It is also noteworthy that the testis is more sensitive to acute stress caused by LPS-induced systemic inflammation. LPS treatment resulted in lower testicular gene expression levels of steroidogenic acute regulatory protein, cholesterol side-chain cleavage enzyme, and cytochrome P450 family 11 subfamily B member 1 after LPS treatment, while no such decrease was found in the adrenal gland. In parallel to the significant decreases in testicular intercellular adhesion molecule 1, tight junction protein 1, and gap junction alpha-1 protein gene expression with LPS treatment, no decrease was found in the epididymis. In the brain, LPS treatment caused higher medial preoptic area (mPOA) activation in the hypothalamus, which is accompanied by elevated blood follicle-stimulating hormone (FSH) and luteinizing hormone (LH) levels, suggesting a disturbed hypothalamic-pituitary-gonad axis function. Besides mPOA, brain c-fos mapping and quantitative analysis demonstrated a broad activation of brain nuclei by LPS, including the anterior cingulate cortex, lateral septum, paraventricular nucleus of the hypothalamus, basolateral amygdala, ventral tegmental area, lateral habenular nucleus, locus coeruleus, Barrington's nucleus, and the nucleus of the solitary tract, accompanied by abnormal animal behavior. Our data showed that LPS-induced inflammation caused not only local testicular damage but also a systemic disturbance at the brain-testis axis level.

Cell-cycle-gated feedback control mediates desensitization to interferon stimulation.

Cells use molecular circuits to interpret and respond to extracellular cues, such as hormones and cytokines, which are often released in a temporally varying fashion. In this study, we combine microfluidics, time-lapse microscopy, and computational modeling to investigate how the type I interferon (IFN)-responsive regulatory network operates in single human cells to process repetitive IFN stimulation. We found that IFN-α pretreatments lead to opposite effects, priming versus desensitization, depending on input durations. These effects are governed by a regulatory network composed of a fast-acting positive feedback loop and a delayed negative feedback loop, mediated by upregulation of ubiquitin-specific peptidase 18 (USP18). We further revealed that USP18 upregulation can only be initiated at the G1/early S phases of cell cycle upon the treatment onset, resulting in heterogeneous and delayed induction kinetics in single cells. This cell cycle gating provides a temporal compartmentalization of feedback loops, enabling duration-dependent desensitization to repetitive stimulations.

Fluorescence Correlation Spectroscopy and Photon Counting Histograms in Finite, Bounded Domains.

Analysis of fluctuations arising as fluorescent particles pass through a focused laser beam has enabled quantitative characterization of a broad range of molecular kinetic processes. Two key mathematical frameworks that have enabled these quantifications are fluorescence correlation spectroscopy (FCS) and photon counting histogram (PCH) analysis. Although these frameworks are effective and accurate when the focused laser beam is well approximated by an infinite Gaussian beam with a waist that is small compared to the size of the region over which the fluorescent particles can diffuse, they cannot be applied to situations in which this region is bounded at the nanoscale. We therefore derived general forms of the FCS and PCH frameworks for bounded systems. The finite-domain form of FCS differs from the classical form in its boundary and initial conditions and requires development of a new Fourier space solution for fitting data. Our finite-domain FCS predicts simulated data accurately and reduces to a previous model for the special case when the system is much larger than the Gaussian beam and can be considered to be infinite. We also derived the PCH form for the bounded systems. Our approach enables estimation of the concentration of diffusing fluorophores within a finite domain for the first time, to our knowledge. The method opens the possibility of quantification of kinetics in several systems for which this has never been possible.